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Cosmic microwave anisotropy satellites promise extremely accurate measures of the 
amplitude of perturbations in the universe. We use a numerical code to test the 
accuracy of existing approximate expressions for the amplitude of perturbations 
produced by single-field inflation models. We find that the second-order Stewart— 
Lyth calculation gives extremely accurate results, typically better than one percent. 
We use our code to carry out an expansion about the general power-law inflation 
solution, providing a fitting function giving results of even higher accuracy. 


1 Motivation 

Two newly approved satellites which will observe microwave background aniso¬ 
tropies, MAP and COBRAS/SAMBA, promise to measure the amplitude of 
these anisotropies at the percent level. Recent theoretical work has demon¬ 
strated that these anisotropies can be predicted from a given mectrum of 
adiabatic perturbations, to an accuracy of one percent or betterld The most 
promising theory for the generation of adiabatic perturbations is cosmological 
inflation, and here we address the question of whether or not it is possible to 
predict the spectrum from an inflationary model at a similar level of accuracy, 
in order to allow one to take full advantage of forthcoming observations. Full 
details can be found in Ref. 2, which also discusses gravitational waves. 

We shall only consider models with a single scalar field (f>, moving in an 
arbitrary potential V{(j)). This is not to say that accurate calculations can¬ 
not be done in more general models, but there typically they have to be done 
on a case-by-case basis, whereas the single-field models can be analyzed si¬ 
multaneously. Our aim is to compute the amplitude Vnik) of the curvature 
perturbation TZ on a given fixed comoving scale k, where the precise terminol¬ 
ogy is defined in Ref. 2. From this one could also compute the spectral index, 
its rate of change, etc. Typically though the corrections we discuss to the 
standard results are observationally negligible except for the amplitude itself. 

2 Framework 

The best calculational framework is that introduced by Mukhanovi and ex¬ 
ploited by Lyth and StewartQ’B The only knowledge we require of the back¬ 
ground (homogeneous) evolution is of the combination z = a(j)/H, where a is 
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the scale factor and H the Hubble parameter. The perturbation can be ex¬ 
pressed by a gauge-invariant potential u = —zTZ. Finally, we use conformal 
time T. Then a fourier mode of the perturbation u obeys the equation 
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Typically we don’t know the background evolution z(t), so we don’t even 
know the equation we are trying to solve. Fortunately though it possesses two 
asymptotic regimes where the solution can be found: 
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The former is the flat-space limit and includes the appropriate quantum nor¬ 
malization. The latter is the growing mode, and corresponds to TZk constant. 
Both regimes are independent of z(t), so the amplitude of perturbations only 
depends on the transition regime. We can therefore expand z(t) about the 
time T where k = aH. 

This expansion_pf the background evolution can be carried out using the 
slow-roll expansion^ Inflation is described using as a fundamental quantity 
the Hubble parameter as a function of the scalar field value, H{(/)). From its 
derivatives one defines a series of slow-roll parameters, the first two being 

The slow-roll approximation demands that all these are small. They can readily 
be computed for a given V {4>). 


3 Results 

We solve the mode equation numerically, in order to compare with existing 
calculations. Figure 1 shows four panels, testing the standard slow-roll approx¬ 
imation r^ult v]l'^ = R^/27r|(^|, the improved second-order result of Stewart 
and I^l|hp an approximation based on using the exact power-law inflation 
result p0 and finally a new expression we derived □ by expanding about a gen¬ 
eral power-law solution. 

We conclude that it is possible to obtain very accurate results numerically, 
and also that anyway the Stewart-Lyth calculation □ is already extremely ac¬ 
curate, being better than two percent accurate in all the parameter region 
currently favoured by observations. 
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Figure 1: A comparison of four analytic expressions with the exact numerical result, with 
the contours showing the relative error. Each point in the e—?? plane corresponds to a 
different inflation model. Above the dashed line is a forbidden region (inflation ends too 
quickly), while the shaded area indicates those models currently preferred by observational 
data constraining the slope of the spectrum and amount of gravitational waves. 
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